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ABSTRACT 

Explaining the observed number counts of submillimetre (submm) galaxies (SMGs) has been 
a longstanding challenge for theoretical models. Surprisingly, recent observations have sug- 
gested that the brightest SMGs are almost exclusively multiple fainter sources blended into 
a single source in the single-dish surveys. This result is in contrast with the predictions of 
our previously presented theoretical model, which includes some effects of blending. In this 
Letter, we consider the implications of an upper limit on the submm flux density for the de- 
mographics of the SMG population. Using a relation amongst submm flux, star formation 
rate (SFR), and dust mass (A/d) from our previous work, we infer the maximum SFR for 
a range of flux density Hmits. For M^ = 2 x 10^ (5 x 10^) Mq, the SFR limit that cor- 
responds to an 870-|J.m flux density (5*870) limit in the range 9 — 12.5 mJy is in the range 
- 630 - 1400 (3600 - 7700) Mq yr'^. The SFR Hmit implies a correspondingly sharp, 
redshift-dependent cutoff in the stellar mass (A/*) function, the value of which we predict 
using the Ss7o-M^, relation predicted by our model. The A/* limit decreases with increasing 
redshift: for an 5*870 limit of 9 - 12.5 mJy, the A/* limit ranges from ~ 4 - 7 x 10^^ M© at 
z^ltO'^3 — 5x 10^^ Mq at z ~ 6. We discuss the few interferometrically detected SMGs 
that may be brighter than the proposed cutoff. Although such objects are certainly interesting, 
inferences based on such objects may not apply to most SMGs. 

Key words: submillimetre: galaxies - galaxies: high-redshift galaxies: luminosity function, 
mass function - galaxies: starburst - infrared: galaxies - cosmology: theory. 



1 INTRODUCTION 

Explaining the observed number counts of high-redshift submil- 
limetre (submm) galaxies (SMGs, which are traditionally de- 
fined as having 850- i i.m flux densit y Sgso S^ 3 — 5 mJy; see 
iBlain et alj 120021 and IShaplevI uOllI for reviews) has long been 
a challenging problem for hierarchical galaxy formation mod- 
els. T hus, SMGs h ave attracted much interest from theoreticians 
(e.g., Baueh et alj 
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Niemi et al.ir2012l ; IShimizu et alJl2012l ; ISomerville et alj|2012l) . To 



explain the observed counts, some have resorted to significant mod- 
ifications of th eir models, such as the use of a top-heavy initial mass 
funct ion (e.g.. lBaugh et alj200^ ; ISwinbank et alj2008l ; lDave et alj 

[2oToh . 

It is becoming increasingly clear that single dish-detected 
SMGs are a heterogeneous population. Although some fraction of 
SMGs show strong evidence for being merger-induced starbursts, 
a significant fraction seem to be early-stage mergers in which the 
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comp onents have projected separation > 10 kpc (e.g.. lEngel et alj 
2010|; for a thorough dis cussion of the observa tional evidence, see 
Havward et al]|2012l) . In lHavward et alj ( 1201 3L hereafter H13), we 
presented a model for the SMG population that includes the fol- 
lowing physically distinct subpopulations: 1. merger-induced star- 
bursts; 2. early-stage mergers in which the galaxies are not yet 
strongly interacting but are blended into a single submm source in 
single-dish surveys, which we term 'galaxy-pair SMGs'; and 3. iso- 
lated discs. Furthermore, there is anecdotal evidence that physically 
and spatially unassociated galaxies blen ded into a single submm 
source also contribute to the population (IWang et alj|201 lb . How- 
ever, the contribution of the last subpopulation is currently very 
poorly constrained, and it is not included in the HI 3 model. 

Except for HI 3, all published theoretical models do not treat 
the effects of blending. In the H13 model, the galaxy-pair SMGs 
contribute significantly to the population at all fluxes, but merger- 
induced starbursts dominate at the bright end. For 5*870 ^ 9 mJy, 
the counts predicted by the H13 model are in good agreement 
with those of single-dish surveys, and the Atacama Large Millime- 
ter/submillimeter Array (ALMA) counts that we pr edicted have 
since been shown to agree very well with observations dKarim et alj 
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l2012h . However, for higher fluxes, the observed ALMA counts de- 
crease significantly more steeply at the bright end than both the 
single-dish counts and the ALMA counts predicted by H13 be- 
cause the brightest (5*870 ^ 9 mJy) single-dish-detected SMGs 
are almost all resolved into two or more sources. This result is 
indeed very surprising. iKarim et al.l (120 12h suggest that their re- 
sults imply an upper limit on t he star formatio n rate (SF R) of an 
SMG of ~ 10^ Mq yr" 1 (for a lSalpeted ( Il955h IMF; for a lKroupal 
( l200lh IMF, the SFR limit is a factor of ~ 2 lower). Based on a 
gas consumption time-scale argument that assumes that SMGs are 
starbursts, they argue that this limit implies that the space density of 
galaxies with gas mass > 5 x 10^" Mq is < 10"'' Mpc"^. In this 
Letter, we elaborate on this proposed limit by 1. using the results 
of dust radiative transfer calculations to infer the SFR limit that 
coiTesponds to a range of flux limits and 2. presenting a redshift- 
dependent stellar mass limit that applies to all constituents of the 
SMG population, not just starbursts. For simplicity, we treat the 
flux density limit suggested by observations as a strict upper limit, 
but in reality, it only implies an upper limit on the number den sity 
of SMGs that are brighter than the limit (see lKarim et al.ll2012h . 

The H13 predictions depend on the assumed stellar mass func- 
tion (SMF), which is especially uncertain at high masses and high 
redshifts. Thus, it is potentially interesting to remove the com- 
ponents of the HI 3 model that depend on the SMF and instead 
work 'backwards', i.e., make inferences about SMG demograph- 
ics from the observations. In this Letter, we use the relationship 
among submm flux density, SFR, and dust mass predicted by our 
hydrodynamical simulations and dust radiative transfer to infer an 
SFR limit from the cutoff flux density suggested by recent ob- 
servations. Furthermore, because our model predicts submm flux 
density as a function of stellar mass and redshift for quiescently 
star-forming galaxies (i.e., isolated discs and the individual com- 
ponents of 'galaxy-pair' SMGs) - and at fixed stellar mass, star- 
burst SMGs should have higher submm flux density than quies- 
cently star-forming SMGs - we can also predict an upper limit on 
the stellar mass of SMGs. 



2 THE STAR FORMATION RATE AND STELLAR MASS 
LIMITS IMPLIED BY THE POTENTIAL FLUX LIMIT 

Dust radiative transfer performed on hydrodynamical simulations 
of isolated disc galaxies and galaxy mergers suggests that the 
submm flux of a galaxy at z ~ 1 — 6 can be parameterized as 
a po wer law of the SFR and dust mass (Md) to within ~ 0.15 
dex l IHavward et al.ll201 iL H13). The physical reason underlying 
the parameterization is that the radiation and dust are in thermal 
equilibrium; the specific form can be motivated by a toy model of 
single-temperature du st that is optically thin to its own emission 
JHavward et al.ll201 ih . The parameterization is relatively insensi- 
tive to the dust geometry and the sub-resolution dust obscuration 
model, and it agrees well with the results of GRASIL JSilva et alj 
1 19981) calculations (A. Benson, private communication). Conse- 
quently, one can infer the SFR of a galaxy from its submm flux, 
albeit with some uncertainty caused by the dependence on Md- 
Solving for the SFR in equations (15) and (16) of H13 (which were 
derived from simulated SMGs with Sgso ^ 15 mJy, so extension 
to higher fluxes is an extrapolation), we have 
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Figure 1. Using the results of our hydrodynamical simulations and dust 
radiative tiansfer, the proposed limit on the submm flux density of SMGs 
IKarim et alJl2012l ; ISmolcic et alj|2012d) can be translated into a limit on 
the SFR. This figure shows the relation between SFR ( Mq yr~^) and 
5'850 (mJy) calculated using Equation (TJ for the M^ values given in the 
legend. The vertical dashed (dotted) line indicates 5*850 = 9 (12.5) mJy. 
The horizontal blue (black) lines indicate the corresponding SFR limits for 
Afd = 2 X 10^ (5 X 10*) Mq, which are ~ 630 (3600) Mq yr'^ for 
5*850 < 9 mJy and - 1400 (7700) Mq yr"! for 5850 < 12.5 mJy. 
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where the subscripts of the SFR variables denote the flux from 
which the SFR is inferred and 5850 and 5i.i are the 850- (j.m and 
1.1-mm flux densities, respectively|jThe relationship between SFR 
and 5850 calculated using Equation lO for various Md values is 
shown in Fig.[T] Because 5850 ~ 2.35i.i (H13), we simply show 
approximate 5i.i values on the top axis of Fig. [T] 

Early interferometric follow-up of SM Gs often focused on 
the bri g htest single dish-detected SMGs (e .g., lYounger et al.ll2007l 
l2008alIbll2009l ; ISmolcic et al.ll201lLl20l2bh . Only recently have in- 
terferometrically observed (sub)mm continuum maps of large (i.e., 
containing more than ~ 10), relatively unbiased samples of less- 
luminou s (i.e., 5850 5, 10 mJy) SMGs become available. Most 
notably. iKarim et alj 1201 21) mapped a sample of 122 SMGs from 
the LABOCA Extended Chandr a Deep Field South Submillimetre 
Survey (LESS ; ' WeiB etal .1120091 ; for the description of the ALMA 
follow-up, see Hodge et al., submitted). At all fluxes, a signifi- 
cant fraction (at least tens of per cent) of the SMGs were resolved 
into two or more components by ALMA. As mentioned above, the 
brightest single-dish SMGs were all resolved into multiple sources; 
none of the resolved sources have Ssto ^, 9 mJy. The i nterfer- 
om etrically imaged sam ples presented in ISmolcic et al.l (l2012al) 
and lBarger et al.l ( 120 12h also tentatively support a steeper cutoff at 
the bright end than expected based on the single-dish counts, al- 
though for these samples, the cutoff occurs at a higher flux density 
(5870 ~ 12.5 mJy), and there are some sources that are potentially 
brighter than the cutoff (see Section|3]for details). The flux density 
at which the cutoff occurs may depend on the field because of cos- 
mic variance (i.e., the highest overdensity sampled can vary from 
field to field) or the different beam sizes of the interferometers used. 
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^ Note that these relations inherently assume a lKroupj 1200 ll) IMF. Also, 
for simplicity, we assume throughout this work that 5860 ^nd 5870 are 
identical to 5850 and use them interchangeably because the conversions 
amongst them are insignificant for our puiposes. 
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Figure 2. The observed sharp cutoff in the SMG number counts suggests 
that there is a correspondingly sharp cutoff in the stellar mass ftinction of 
star-forming galaxies at a given redshift. The above plot shows the predicted 
Mir{z) limit that corresponds to an 5g7o hmit in the range 9 — 12.5 mJy, 
where a higher 5870 hmit corresponds to a higher A/* limit. See the text 
for details regarding the calculation of the limit. 



^Because submm flux depends weakly on SFR JHavward at alj 

120111) . the inferred SFR limit is sensitive to the maximum value 
of the submm flux observed. The SFR limit is also sensitive to 
the value of Md. The Md values inferred from fitting the far- 
infrared spectral energy distributi ons of SMGs can be i n excess of 
10^ M0 but are often lower (e.g. jMagneni et al.ll2012h . However, 
the accuracy to which Md can be determined is still debated (e.g., 
lOale et al.ll2012h . Thus, we will determine the SFR range that cor- 
responds to 5870 cutoffs in the range 9 — 12.5 mJy for a range of 
Md values that are c haracteristic of those inferred for SMGs (e.g., 
lMagnemetaLll2012h . Fig. [U indicates the Ssro = 9 (12.5) mJy 
limit with a vertical dashed (dotted) line. The corresponding SFR 
limits for Md = 2 x lO'-* (5 x 10*) M© are shown as blue (black) 
horizontal lines. Assuming Md ^ 2 x 10^ (5 x 10**) Mq, Ss7o 
limits in the range 9 — 12.5 mJy correspond to maximum SFRs in 
the range ~ 630 - 1400 (3600 - 7700) M© yr~\ Note that the 
SFR limit for Md = 2 x lO'' M© is similar to that proposed by 
lKarimetalJJ2012h . 



2.1 Quiescently star-forming SMGs 

In the H13 model, both isolated disc galaxies (which dominate 
the faint, i.e., Ss50 < 2 mJy, SMG population) and the galaxy- 
pair SMGs (i.e., early-stage mergers in which the discs are not yet 
strongly interacting) form stars quiescently (as opposed to via the 
starburst mode). The H13 relations for the redshift evolution of the 
gas content, size, and metallicity of galaxies can be used to predict 
a quiescently star-forming galaxy's SFR and Aid as a function of 
its M* and z. Because the (sub)mm flux of our simulated galaxies 
can be parameterized as a function of SFR and Md alone, the H13 
model can be used to determine the (sub)mm flux as a function 
of Af, and z (see equations 17 and 18 of H13). Thus for a given 
redshift, the SMG selection effectively selects the most massive 
quiescently star-forming galaxies on the main sequence because 
the Ss5o-M* relation predicted by our model for fixed redshift 
is monotonic (although in practice, scatter in the SFR-M<, rela- 
tion w ill cause scatter in the Sgso-M* relation; |Michalowski_et_alJ 
l2012h . Therefore, the observed flux density limits can be translated 
into upper limits on M* as a function of z. 



Fig. m shows the inferred Af, limit versus z for an 5870 limit 
in the range 9 — 12.5 mJy. This limit is calculated by combining 
the redshift-dependent SFR-Af* relation predicted by our model, 
which is tied to the dependence of galaxy gas fractions and sizes 
on A/* and z, and the dependence of the dust mass on A/* and z, 
which depends on the gas fraction and metallicity evolution. Then, 
the flux limit and the A/* and A/d functions are substituted into Eq. 
(TJ, and the resulting equation is solved for Af*. Because our model 
predicts an SFR-A/* relation that shifts upward as z increases, the 
mass limit decreases with redshift. Thus, it is natural to expect that 
the brightest quiescently star-forming SMGs are at higher redshifts 
than the median for the SMG population if the high-mass end of the 
stellar mass function evolves less rapidly than the mass limit plot- 
ted in Fig. [21 n Indeed, we shall see below that most of the very 
few sources with 5870 ^ 12.5 mJy have redshifts significantly 
greater than the median redshift of the population. The M* limit 
that coiTesponds to an 5870 limit of 9 — 12.5 mJy ranges from 
~ 4 - 7 X 10^^ Mq at 2 ~ 1 to ~ 3 - 5 X lO" M© at z ~ 6. 

It is also interesting to consider the specific SFR, SSFR = 
SFR/A/*. In the H13 model, the M* and z dependence of the SSFR 
are determined by the dependence of galaxy gas fractions and sizes 
on M* and z. At fixed redshift, SSFR decreases with M^. ; the scal- 
ing is similar to that of the star-forming galaxies in the lKarim et alj 
( I2OIII) sample, SSFR ex M~°-^. At fixed M^, SSFR increases with 
redshift. For example, a galaxy with A/* — 10^^ M©, approxi- 
mately the minimum Af* necessary for a main sequence galaxy to 
contribute to the SMG population (see fig. 3 of HI3), has SSFR 
that ranges from 0.4 Gyr^^ at 2 ~ 1 to 6.7 Gyr^^ at 2 ~ 6. 
Consequently, we expect the higher-redshift main sequence SMGs 
to have higher SSFRs. Furthermore, the H13 model suggests that 
main sequence SMGs should have SSFR < 7 Gyr~^. 



2.2 Merger-induced starburst SMGs 

In the H13 model, merger-induced starbursts contribute signifi- 
cantly to the SMG population, and they dominate the bright (i.e.. 
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5 mJy) end of the population. Merger-induced starbursts 



can be outliers from the main sequence (i.e., have higher SFR for 
a fixed Mi,, but not all merger-induced starbursts are outliers from 
the main sequence ;|Hay ward et al. 2012). Indeed, observations in- 
dicate that SMGs are a mi x of main sequence members and outliers 
JMichalowski et al.ll2012r) . (Note that, as discussed below, some of 
the SMGs that are outliers from the main sequence could also be 
gravitationally unstable discs.) If the sources with flux density near 
the proposed cutoff flux density all lie significantly above the SFR- 
A/* relation, then the upper limit on M, presented in Section 12.11 
would overestimate the true M* limit. 

In the merger simulations of HI 3, which can account for 
the observed abundance of SMGs that are fainter than the pos- 
sible flux limit, the highest SFR attained is ~ 5400 M© yr^^. 
The peak SFR is typically limited by gas consumption rather than 
AGN feedbaclij; even without AGN feedback, an SFR in excess 



■^ This conclusion will not hold at redshifts beyond which the nor- 
malisation of the SFR-Af* relation no longer increases. Observation- 
ally, the normalisation of the S FR-M<, relation seems to increase wit h 
redshift out to 2 ~ 2 (e.g., iNoeske et al] |2007| ; iDaddi et all l2007t) ; 
whether t he trend continues to higher redshift is s t ill an open ques- 
tion (e.g.jGonzalez et al]|2010l;lRodighiero et aljboiol : iKarim et alj|201 ll ; 
iMcLure et al.l201 iMWhitaker et al.l20l'2l\ 

■^ Note that the H13 simulations do not include starburst -driven winds. 
However, extremely high-resolution merger simulations ( IHopkins et alj 
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Table 1. Sources that potentially exceed the flux density limit 



Identifier{s) 



Flux density" (mjy) Reference* 



GN20 22.9 ± 2.8 (1) 

AzTECS, Cosla-73, AzTEC/C2 21.6 ± 2.3 (2) 

BR1202-0725 North 18.8 ± 2.8 (3) 

BR1202-0725 South 18.0 ± 2.7 (3) 

MMl,Cosbo-l 16.8 ±1.5 (4) 

AzTEC4, AzTEC/C4 14.4 ± 1.9 (5) 

AzTECl,COSLA-60, AzTEC/C5 13.8 ± 2.3 (5), (6) 

LH 850.02 13.4 ± 2.4 (7) 

AzTEC2, COSLA-4, AzTEC/C3 12.4 ± 1.0 (5) 

AzTEC12, AzTEC/C18 12.8 ± 2.9 (8) 

AzTEC/Cl,COSLA-89 12.4 ± 3.7 (9) 

AzTEC7 12.0 ± 1.5 (5) 

GNIO, GOODS 850-5 12.0 ± 1.4 (10) 

MM J 120546-074 1.5 6.5 ± 0.8 (11) 
"For all but MM J120546-0741.5, the flux density is for a wavelength in 
the range 850 — 8 90 |J-m; for that s our ce, the 1.2-mm flux d ens ity is given. 

■ Referenc e(s): (l)llono et al.ll200d (2)lYounger et al.ll2009l (3) IWagg et al ' 



201 d (4) lAi-avena et al J |20 1 Ol (5) lYounger et al.[ 20071 (6) TYounger et al. 
2008b[. (7) lYounger e tal.ll2008j. ( 8) lYou nger et al." 2009l (9 ) ISmolcic etal 
2012bl. (10) Barger et al.l2012l (ll) lDann erbauer et al.l2002i 



of 10'^ Mq yr^^ can only be sustained for a few x 10 Myr unless 
fresh gas is supplied to the galaxy at a comparable rate. Thus, the 
potential SFR limit discussed in this work m ay simply be a con - 
sequence of limited gas supply, as argued by JKarim et alj j2012h . 
rather than evidence for strong feedback during the starburst. 

Also, note that the starburst driven at coalescence in major 
mergers can elevate the SSFR beyond the aforementioned maxi- 
mum for quiescently star-forming galaxies. In simulations of merg- 
ers with baryonic mass ratio > 1/3, the SSFR can be boosted by 
~ 10 — 20 X for a very short period (~ 30 Myr at most) and a factor 
of a few for as long as ~ 100 Myr. 



3 SOURCES THAT POTENTIALLY EXCEED THE FLUX 
DENSITY LIMIT 

We have performed a literature search to identify SMGs for which 
interferometric continuum imaging indicates Sgro ii 12.5 mJy; 
they are especially interesting targets for follow-up observations to 
determine whether they are indeed brighter than the proposed flux 
limit and to discern the physical nature of these potentially extreme 
sources. We are aware of seven sources at 2 > 1 (we are uninter- 
ested in low-z interlopers) that are at least \a and eight additional 
sources that are possibly brighter than the 5*870 ~ 12.5 mJy cutoff. 
The names of these sources, their flux densities, and references to 
the interferometric continuum observations of which we are aware 
that confirm the high fluxes are presented in Table [T] 

Many of the sources that are potentially brighter than the 
proposed flux limit are from the AzTEC survey of a subfield 
of the Cosmic Ev olution Survey (CO SMOS; IScoville et al.ll2007h 
fiel d performed by IScott et al.l ( 120081) . The number counts inferred 
by Aretxaga et al.l OOllI) for a larger subfield that contains the 
IScott et al.l ( 120081) field are significantly higher at the bright end 



I2OI2I) that include sophisticated presciiptions for feedback from super- 
novae, stellar winds, radiation pressure, and photoionisation suggest that 
the peak SFR of the merger-induced starbu rst is relatively unaffected by 
starburst-driven winds JHopkins et al.l2013bl) . 



than those derived from the similar-sized Submillimetre Com- 
mon User Bolometer Arra y (SCUBA) Half Degre e Extragalac- 
tic Survey (SHADES) field ( JAustermann et alj|2010l) . The brighter 
sources tend to coincide with overdensities of galaxies at z < 1.1, 
which may suggest that ga laxy-galaxy lensing is responsible for the 
exces s of bright sources ( lAustermann et al.ll2009t JAretxaga et alj 



l201in . Some of the other brightest sources may also be lensed. 

However, it is unlikely that lensing is responsible for all of the 
sources that are potentially brighter than the proposed flux density 
limit, and it is worthwhile to consider whether the sources that are 
brighter than the proposed limit are exceptional in some manner. 
As we argued above, many of the brightest sources are at redshifts 
significantly greater than the median redshi ft of the SMG popu- 
lation (^ 2.6 - 3.5; e.g., Yun et al. 2012i; ISmolci£"e t al. 2012^; 
IWeiB et alj|2013l) . For examp le, GN20 and GNI O are both located 
in a z = 4.04 protocluster JPaddi et al.ll2009allbl) , AzTECl is at 
z = 4.6 dSmolcic et ai]|201lh . and BR1202-0725 North and South 
are in a z = 4.69 overdensity. The scarcity of bright sources at 
higher redshifts may be c aused by a lack of sufficiently massive 
galaxies at those redshifts ( lHavwardll2012l) and/or insufficient dust 
content because of the li mited time available to produce dust (e.g., 
JMichalowski et alj20ld) . If indeed the bulk of the sources brighter 
than the proposed cutoff are at z > 4, this result may suggest that 
the normalisation of the SFR-M* relation continues to increase to 
at least that redshift; otherwise, there would be no clear reason for 
the brightest sources to be preferentially at z > 4, and the evolution 
of both the stellar mass function and dust content of galaxies should 
tend to make the brightest SMGs less abundant at those redshifts. 

As an example of what can potentially be learned from de- 
tailed follow-up of the sources that are potentiall y brighter than th e 
proposed cutoff, we consider the case of GN20 ( jPope et al.ll200q) . 
which is an especially interesting and well-studied source. Re- 
cent high-resolution interferometric observations performed with 
the Very Large Array (VLA) indicate that GN20 is a very extended 
(~ 14 kpc in diameter) gas-rich disc that contains at least five 
molecular gas clumps, each of which contain a few percent of the 
total molecular gas mass of the galaxy (Hodge et al. 2012). In our 
model, the redshift evolution of the SFR- A/* relation alone cannot 
account for the extrem ely high flux densi ty of GN20 (which does 
not seem to be lensed; ICarilli et aLJuOlOl) : Fig. |2] indicates that at 
z -- 4, a galaxy with M* < 5 x lO" M© would have Sgro < 9 
mJy, and the inferr ed M* value for GN20 (M* ~ 2 x lO" Mq; 
JDaddi et alj l2009bl) is less than this value. Thus, GN20 seems 
to be an outlier from the main sequence (see als o JDaddi et alj 
l2009bl ; ICarilli et alj|2010l I2OIII; JHodge et alj|2012l) : its SSFR is 
~ 8.6 Gyr~^ dMagdis et alj2011r) . but in our model, a z ~ 4 main 
sequence galaxy with M^, = 2x lO" Mq has SSFR = 2.3 Gyr"\ 
~ 4x less than that of GN20. The simplest interpretation of GN20 
is that it is an extremely massive, gas-rich, gravitationally unsta- 
ble disc undergoing rapid fragmentation, but whether such galaxies 
constitute a significant fraction of the SMG population is unclear. 



4 SUMMARY AND DISCUSSION 

Recent interferometric imaging of single-dish-detected SMGs has 
demonstrated that most of the brightest (Ssro ^ 9 — 12.5 mJy) 
sources are actually multiple distinct sources blended into a single 
source. In this Letter, we have used the relationship among submm 
flux density, SFR, and dust mass predicted by our combination of 
hydrodynamical simulations and dust radiative transfer and the red- 
shift evolution of the SFR-M* and Md-M* relations predicted by 
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the H13 model to predict limits on the SFRs and stellar masses of 
SMGs from the sharp cutoff in the number counts suggested by re- 
cent observations. For Md = 2 x 10® (5 x 10*) Mq, an Sgro limit 
in the range 9 — 12.5 mjy corresponds to an SFR limit in the range 
~ 630 - 1400 (3600 - 7700) M© yr" ^ The Af * limit ranges from 
~ 4 - 7 X 10^^ M© at 2 ~ 1 to ~ 3 - 5 X lO" Mq at z ~ 6. 

Follow-up observations of the sources that are potentially 
brighter than the proposed flux limit are crucial for understanding 
whether these sources are exceptional in some sense (besides their 
bright submm fluxes). Are they predominantly merger-induced 
starbursts, gravitationally unstable discs, or something else? Such 
observations may yield new insight into star formation and feed- 
back at high redshift. 

However, one should be cautious when attempting to make 
inferences about the SMG population as a whole from observa- 
tions of potentially extreme objects such as GN20, which is one 
of the brightest SMGs known and, contrary to naive expectations, 
is an extended very gas-rich disc galaxy rather than a starburst 
([Hodge et al. 20121). Because the H13 model, which does not in- 
clude violently unstable discs, can largely account for the ob- 
served SMG number counts, it is not clear that the SMG population 
must include a significant number of such objects. Even if some 
SMGs are massive, gas-rich, gravitationally unstable discs, this 
does not directly imply that they are fueled by ' cold-mode accre- 
tion' JBimboim & Dekell2003r Keres et alj2005l), as has been sug- 



gested by some au thors (e.g.. lDave et al.ll2010tlCarilli et alJuOld : 
iHodge et alj|2012r) . Rather, it simply requires that the discs are suf- 
ficiently gas-rich to be gravitationall y unstable and is agno stic re- 
garding the method of gas supply (see iHopkins et al.l2013al for fur- 
ther discussion). 

Finally, we note that we have not attempted to reconcile the 
discrepancy betwe en the counts predicted by H13 and those of 
iKarim et alj ( 120121) in this work. One potential reason for the dis- 
crepancy is that in HI 3, we did not include the effects of blending 
of physically and spatially unassociated galaxies. We also did not 
include mergers of more than two galaxies. These potential contri- 
butions will be modeled in future work. 
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